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Abstract 

A consistent procedure for regularization of divergences and for the subsequent renor- 
malization of the string tension is proposed in the framework of the one-loop calculation 
of the interquark potential generated by the Polyakov-Kleinert string. In this way, a jus- 
tification of the formal treatment of divergences by analytic continuation of the Riemann 
and Epstein-Hurwitz zeta functions is given. A spectral representation for the renormal- 
ized string energy at zero temperature is derived, which enables one to find the Casimir 
energy in this string model at nonzero temperature very easy. 
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1 Introduction 



A consistent method to treat the divergences in quantum field theory is known to be the 
following [1]. Divergent expressions must at first be regularized, for example, by the Pauli- 
Willars method, then subtractions justified by transition to physical (observable) parameters 
of the theory should be done. After that the regularization is to be removed. 

Together with this approach there are widely used methods that do not apply explicit 
regularization and renormalization but which nevertheless give finite answer. First of all, it is 
the zeta function technique. The main idea of this approach is the following [2-4]. One assumes 
that the divergent sum J2 n u n of eigenvalues of the operator^ determining the dynamics in the 
model under consideration is equal to the value of the zeta function for this operator, C( s )? when 
s — > —1. At first the function £(s) is defined by the formula ((s) = I]„c<;~ s for Res > 1, and 
then it is analytically continued to Re s < 1 possibly save for isolated points. In the case of the 
Dirichlet boundary conditions for the two-dimensional Laplace and Helmholtz operators this 
function appears to be the Riemann ^-function or the Epstein-Hurwitz ^-function, respectively. 
These functions are widely used in calculations of the Casimir energy in field [4,5] and string 
models [6]. 

Undoubtedly such a formal method to treat divergences needs justification in each particular 
case [7,8]. The more so, there are examples when analytic continuation leads to ambiguities [9]. 
To justify this approach, it is necessary to show that it gives the same results as the standard 
renormalization procedure with regularization and subtraction. It is this problem that will 
be considered in the present paper in the framework of one-loop calculation of the interquark 
potential (or the Casimir energy) in the rigid string model. This model is chosen because here 
both the Riemann and Epstein-Hurwitz ^-functions are employed. 

The interquark potential generated by a rigid string was studied in a number of papers 
by making use of the perturbation theory and variational estimation of the functional integral 
(see, for example, Ref. [10] and papers cited therein). These results are well-known. Therefore 
attention will be basically paid to development of the consistent procedure of renormalization 
and to justification, on this basis, the results obtained by ^-function method. 

The layout of the paper is as follows. In Section 2, the interquark potential generated by 
a rigid string is calculated in the one-loop approximation, the standard method of analytic 
continuation of the Riemann and Epstein-Hurwitz ^-functions being used. In Section 3, the 
consistent regularization of the divergences and the string tension renormalization are carried 
out. Unlike the ^-function method, the finite expression for the string potential is derived 
here uniquely. Moreover in our approach the renormalized string energy at zero temperature is 
obtained in terms of the spectral representation that can be directly generalized to a finite tem- 
perature. In the Conclusion (Section 4), the obtained results are discussed in short. Auxiliary 
material concerning the details of the calculations is given in Appendices A and B. 



^^Most commonly this operator is the Laplace operator (—A). The zeta function regularization is usually 
applied to the Euclidean version of the models where one has to do with elliptic operators [3] . 
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2 Interquark potential generated by rigid string in one- 
loop approximation 

We consider the most simple example of the application of the Riemann and Epstein-Hurwitz 
^-functions. This is the calculation of the interquark potential generated by the Polyakov- 
Kleinert string [11,12] in the one-loop approximation [13]. In spite of its simplicity, this example 
demonstrates the main features of the approach. 

For our purpose the quadratic approximation to the Polyakov-Kleinert string action is suf- 
ficient 

P R r , / v -, 

S" = M*fdtfdT l + -u 1-^A (-A)u . (2.1) 
o L V o / 

Here Mg is the string tension, u(t, r) = (w 1 (t, r), u 2 (t, r), . . . , u D ~ 2 (t, r)) are the transverse string 
coordinates in D-dimensional space-time, and a is a dimensionless parameter characterizing 
the string rigidity, a > 0, R is the distance between quarks connected by string, i.e., the 
string length. The Euclidean action is considered, therefore the operator A in ( |2.1| ) is the 
two-dimensional Laplace operator A = d 2 /dt 2 + d 2 /dr 2 . The "time" variable t ranges in the 
interval < t < f3, where (3 = 1/T is the inverse temperature. 

The action (|2.1|) should be completed by boundary conditions for string coordinates at the 
points r = and r = R. Usually a string with fixed ends is considered 

u(t, 0) = u(t, R) = . (2.2) 

This corresponds to the static interquark potential. The string potential V(R) is defined in a 
standard way 

exp[-/3V(R)] = J [Du] exp {-^[u]}, j3 -»• oo. (2.3) 

The functional integral in ( ^.3| ) is taken over string coordinates u(i, r) that satisfy periodic 
conditions in the time variable t 

u(t,r) = u(t + P,r). (2.4) 
Inserting (|2.1|) into ( p.3|) one obtains after the functional integration when j3 — > oo 

V(R) = M 2 R + ^^Tr ln(-A) + ^^Tr In (1 - A j . (2.5) 



2p y ' 2(3 I M ( 



For calculating the traces in ( |2.5p the eigenvalues of the operators (—A) and [1 — («/Mq)A] 
with the boundary conditions (|2.2|) and the periodicity conditions ( |2.4j ) are needed 

( A)(^ nm ^nm^nmi 

(1 - -^A j = (2-6) 



M 2 ' 



Using the Fourier expansion we find [13] 



Anm — + ^ 2 



3 



U = + (2-7) 

where Q„ = 2irn/f3, n = 0, ±1, ±2, . . . are the Matsubara frequencies, u m = mir/R, m = 
1,2,... are positive roots of the equation 

sm(uR) = 0, (2.8) 



and Uk = v {kTr/R) 2 + M^/a, k = 1, 2, . . . are those of the equation 



sin [R\JCj 2 - M§/aj = 0. (2.9) 

Summation over the Matsubara frequencies Q n can be accomplished by making use of the 
known methods [13]. Upon taking the limit (3 — ► oo, the potential generated by the string 
assumes the form 

V(R) = MqR +(D- 2)(E { o ) + E { c ] ), (2.10) 

where E& are the Casimir energies corresponding to both the modes of the rigid string 

oscillations 



E$ = — Y n (2.H) 
c 2R 2^ I ) 



oo 



Summation of the divergent series ( 2.1 lj ) and (p. 12 ) by analytic continuation of the ^-function 



is now commonly used. Nevertheless we remind the main steps of this approach in short. 
We begin with the first sum Q2.ll ). According to the scheme outlined in the Introduction, 



we first have to consider the function 

C(s) = E™" S > Res>l, (2.13) 



oo 



71=1 



and then to continue it analytically to the region Re s < 1. In this case ((s) is the Riemann 
^-function. Analytic continuation of the function (|2.13| ) to the rest of the complex plane s, 
with the exception of the point s — 1, is performed by the contour integral [14] 

c 

where the contour C is shown in Fig. 1. This contour should avoid the points z = ±2nni (n = 
1, 2, 3, . . .). Because of the multiplier T(l — s) in Q2.14 ) the Riemann (^-function has a simple 
pole at s = 1 with the residue equal to 1 

C( s ) = ^- + 7 + ..., (2.15) 
s — 1 

where 7 is the Euler constant 

N 



7 

N 



lim ( Y --Iiin] . (2.16) 
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The (^-function defined by integral (|2.14|) satisfies the reflection formula [14] 



C(l -s) = 2 (2n)-° cos (y) r(s)((s). (2.17) 

According to the scheme outlined above we have to attribute the value £(— 1) to the sum of 
the divergent series (|2.11 ) . For s = — 1 the integral representation ( 2.14Q gives 



c 

Since the integrand is single-valued in the plane z, the integration contour in Fig. 1 can be 
closed. As a result, ((— 1) is equal to the integrand residue at z = 0. To find the residue, we 
can use the definition of Bernoulli numbers [14] 

* I t 2 t A 

l--t + B 1 --B 2 - + ..., 



e'-l 2 x 2 z 4! 

where B 1 — 1/6, B 2 — 1/30, — Thus we obtain 

C(-D = -fti = -4 (2.19) 

Finally we attribute the value 

4 1 ) = _L V n = — C(-l) = -— . (2.20) 
c 2i2^ 2iT v ; 24i? v ; 

to the sum of the divergent series ( |2.11| ). In the theory of divergent series [15] this summation 
of the series ( |2 . 1 1| ) is referred to as the Ramanujan summation. Obviously, this method is not 
universal. For example, it cannot be applied directly to a divergent series Y,™=i n ~ 1 because 
the zeta-function has a pole at the point s = 1 (see Eq. (|2.15| )). Admitting the convention 
about the rejection of the pole singularity, as it is usually done in the analytic regularization 
met ho cfl, we get 

From ( |2.15|) , (|2.16|) it follows that the pole of the Riemann (^-function at the point s — 1 is 
responsible for the logarithmic divergences. 

Summarizing, we arrive at the conclusion that the Riemann ^-function method enables 
one to obtain the finite value of the Casimir energy ( |2.11 ) without explicit regularization, pole 
singularity rejection, and explicit renormalization. However in the case of divergent series (|2.12|) 
the analytic continuation technique requires additional assumptions. 

To sum the series ( |2.12|) , we have to consider the Epstein- Hurwitz zeta-function (eh(s,p) 
defined by the formula 

OO 

( E H(s,a 2 ) = J2(n 2 + a 2 )- s , (2.21) 
n.=i 



2 It is implicitly assumed that the singularities are taken away by renormalization of parameters in the theory 
under consideration. 
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where s > 1/2. Let us remind briefly how to accomplish an analytic continuation of this series 
to the region s < 1/2. 

Using the integral representation for the Euler gamma function [14] 



(n 2 + a 2 )- s T(s) = J t^e-^+^dt, 
o 

we can replace each term in series ( |2.21| ) by the integral 

(EH(s,a 2 ) = -±- lt^±e-^) d t. 

1 \ S ) J Q n =l 

The Jacobi ^-function appearing in ( |2.23| ) 9(t) = J2™=i e - " 2 ' has the property [14] 
Substituting JT241 ) into (^23|) we obtain 

^r( s -i/2) 



(2.22) 



(2.23) 



(2.24) 



Ceh(s, a 



r s) 



n=l 



+ 



s-3/2 



2T(s) 



2^2' 



tc n 



dt. 



(2.25) 



The multiplier exp(— ta 2 — n 2 n 2 /t) ensures the convergence of the integral in ( |2.25|) for all 
s. Then this integral is expressed in terms of the modified Bessel function K v (z) having the 
integral representation [16] 



J x v ~ l exp (-1 - 5x\dx = 2 (?) K v {2^8) 



(2.26) 



KJz) = KJz). 



In the case under consideration, 7 and 8 are positive quantities: 7 = ir 2 n 2 , 5 = a 2 . Now we can 
rewrite formula (|2.25|) as follows 



,2\-s 



(eh(s, a 
2tt s 



^(S - 1/2) 2N- S +l/2 + 



+ 



r s) 



2 2r(s) 

,/2+l/4g n ,-l/2^_ i/2(27rnv ^ 



(2.27) 



n=l 



The series obtained converges for all s as the modified Bessel function has the asymptotics [16] 



1/2 



z — >■ 00. 



(2.28) 
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Therefore, the singularities of the function (eh(s, a 2 ) are due to the singularities of T(s — 1/2) 
in (2.27), i.e., (eh(s, a 2 ) has first order poles at the points 



.= 5. -5. (2.29) 

Thus formula ( ^.27[ ) affords an analytic continuation of fl2.21| ) to the region s < 1/2 except for 
the points Q2.29|) . Since (eh(s, « 2 ) has a pole at s = —1/2, function (|2.27| ) can be used for 
obtaining only the regularized Casimir energy E$ Teg 

Ef „, = _M» _ MSR _ M* g n ., (2nM^R\ 

In order for the regularization to be removed the first and second terms in the right-hand side 
of ( ^.3U[ ) should be omitted (see, for example, [17-20]) 



p (2)ren M ^ _ 2nM Q R \ 

E c =~^L n K i\ 7^- I ■ ( 2 - 31 2 

2 V« n=l V V« 



Rejection of the second term in (|2.30|) proportional to T(— 1) is natural in the analytic contin- 
uation methodP]. As for the first term — Mo/(4y/a), its rejection seems to be rather arbitrary. 
Usually this is motivated by the fact that this term is independent of R and, as a consequence, 
does not contribute to the Casimir force. However, this argument does not explain the rejection 

( 2") 

of the i?-independent term in the interquark potential, i.e., in E c (R). In the general case those 
terms may be essential for the description of quark-quark interaction inside hadrons. Only the 
consistent renormalization with preliminary regularization and subsequent subtraction can jus- 
tify the rejection of both the first and second terms in the right-hand side of ( |2.30| ). This will 
be demonstrated in the next Section. 



3 Renormalization of the string tension and removal of 
the divergences 

Let us calculate the interquark potential (|2.5|), Q2.10| ) applying the standard renormalization 
technique. The initial model includes two parameters: string tension Mq and a dimensionless 
constant a characterizing the string rigidity. In the one-loop approximation, only the string 
tension is renormalized. 

The renormalized potential of the string at large distances should coincide with its classical 
expressionF] 

V ren (R) l^oe = M 2 R, (3.1) 

3 If fields are considered in a bounded space region, then this procedure is interpreted as the subtraction of 
infinite space contribution [21]. 

4 In the framework of a string model, the potential linearly rising at large distances is the classical string 
energy considered as a function of its length R when R — > 00 [6]. On the microscopic level (QCD level), the very 
appearance of the collective string degrees of freedom is interpreted as a complicated nonperturbative effect in 
quantum dynamics of gluon and quark fields closely related to nontrivial properties of the QCD vacuum. 
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where M 2 is the renormalized string tension, whose explicit expression for be obtained fur- 



ther. Starting with (|2.10|) and taking into account the necessity to regularize all the divergent 
expressions, we represent V ren (R) as 



V ren (R) = M 2 R+(D-2) 
= M 2 R+(D-2){ 



E$ re9 (R,A) + E ( c z>re9 {R,A) 



i(2) reg. 



A^oo 



E$ re9 (R,A) + E f 



g )re9 (R,A) 



E$ )re9 (R 



00, A) + E c 



+(D - 2) 



E$ re9 (R 



M 2 R + (D - 2) E§ ) ren {R) + E%> ren (R) 



(2)re9 (ii!^oo,A)]} 
oo,A) + 4 2) re9 (R 

-i(2) ren , 



A^oo 

00, A) 



+ 



(3.2) 



where A is a regularization parameter; M 2 is the renormalized value of the string tension, 



M 



M 2 + 



D 



R 



E, 



(1) reg 
C 



(R -> 00, A) + E { c ] re9 {R -> 00, A) 



A— >oo 



and E, 



(i) ren 

C ■> 



1, 2 are the renormalized Casimir energies fl2.11|) and (|2.12 



E { § ren (R) 



E { i ] re9 {R, A) - E%> re9 (R -> 00, A) 



(i) reg/ 



A^oo 



1, 2. 



(3.3) 



(3.4) 



To regularize divergent series ( |2.11| ) and ( |2.12j ), we substitute them by finite sums that can 
be represented in terms of the Cauchy integrals [14] 



- -' I /(*) 



k I 



(3.5) 



Here f(z) is an analytic function having, in a region surrounded by contour C, zeroes of order 
rik at points z = and poles of order pi at points z = b\. As a function f{z) we substitute the 
right-hand sides of frequency equations (|2.8|) and (|2.9|) into ( |3.5|) and choose the contour C so 
as to include first positive roots of the corresponding equations. Functions ( ^.8|) and ( |2.9| ) 
have zeroes of the first order on the real axis and have no poles. Therefore only the first sum 
with rik = 1 remains in the right-hand side of Q3.5Q . 

First, we obtain the regularized Casimir energy ( |2.11| ) 



E, 



(1) reg 
C 



Am J sin (uR) 

c v ' 



(3.6) 



where the contour C is shown in Fig. 2. All the singularities of the integrand in (3.6) being 
situated on the real axis, it is possible to deform the contour C to C continuously (see Fig. 2) 
Now the regularization parameter is the radius A of the semicircle entering into the contour C . 

To determine the counterterms according to (|3.1| )- (|3.4j ), it is necessary to find the asymp- 
totics of E^ re9 (R) for R — > 00 and fixed A. On the semicircle of radius A (Fig. 2) the 
asymptotics of the integrand for R — ► 00 is the integrand itself because of its oscillating char- 
acter. Consequently, the result of integration along this part of the counter C is completely 
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absorbed by the counterterm and does not give any finite contribution to E^ ren (R) . Now let 
us turn to the integral along the interval (— iA, iA) on the imaginary axis 



• i l.i <•■■■/, •„ , R f ... cosh(%) 



E^ {R ,K ) = -- jy -^ )dv . (3.7) 



To find the asymptotics needed, we integrate in (3.7) by parts 

A A 

E§ )re9 {R,A) = / y d(\n\smh(Ry)\) = ln[mnh(Ai*)] + — / dy\n[sinh(Ry)]. 

An J 2tt 27r J 



When R — > oo 



A 

E$ reg {R -> oo, A) = -— ln[sinh(Ai?)] + ^- f{Ry- In 2) dy. (3.8) 

27T 27T J 



Inserting (3.8) into (3.4) we obtain the finite value for the renormalized Casimir energy (2.11) 

oo oo 

4" -(*) = ± J in (1 - «L, = -f / (3.9) 



The last formula is derived by integrating by parts. It is interesting to note that (3.9) is 
expressed in terms of the value of the Riemann ^-function at the point s = 2. Really, 



DC 



r uoduo r Lue~ UUJ , ^ r nrut , L(2) ^ 1 1 . , 

/ 7 = / du = Y j uje~ anu} duj = -V E -? = tC(2). 3.10 

i e"-l 7 1 - e" aa; 4^ 7 a 2 n 2 a 2 w v 7 

n_i n ~ i 

In view of this, Eq. fl3.8|) can be rewritten as 



E g)- (jR) = __i_ c(2) _ (3n) 

Thus, under consistent renormalization, the sum of divergent series ( |2.11 ) is also defined through 
the Riemann (^-function, but now another range of its definition is used, namely, the region 
Res > 1. Here ((s) is defined by convergent series ( |2.13| ). 

With the help of the Riemann reflection formula ( 2.17 ) the value of ^-function at s = 2 
entering into (3.11) can be expressed through ((— 1) 

C(2) = -2vr 2 C(-l). (3.12) 

Final renormalized formula for the Casimir energy Q2.11D assumes the same form as that 
obtained by analytic continuation of the Riemann ^-function ( Eq. ( |2.2U[ ) ) 

(1) ren {R) = JL«-1) = — ( -—) = — . (3.13) 

Thus, there is a complete agreement between two outlined approaches to the calculation of 
the finite value of . 
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Before we turn to consideration of the series ( 2.12|) , let us make a short remark concerning 
formula ( [3.9|) . Discarding the minus sign, the integrand in ( |3.9|) has a form of the Planck energy 
distribution in the spectrum of one-dimensional black-body with temperature 1/2R. 

The renormalized value of the Casimir energy E§ ] (see Eqs. (gig) , (|J)) can be obtained 
in the same way as it was done above. Substitution of the frequency equation ( |2.9| ) into (3.5) 
gives 

,2 



(2) reg 
C 



R r cos( 



c 



LUn 



10 



sin 



or 



do;, 



(3.14) 



where Uq = M^/a. When choosing the contour C one has to take into account the branch 
points of the integrand u = ±c<Jo. To select the single-valued branch of the function, we 
make a cut connecting the branch points along the real axis. After that the contour can be 
chosen as shown in Fig. 3. Again integration along the semicircle of radius A contributes only 
to the counterterm. The integrals along the edges of the cut are mutually cancelled, and the 
contribution I\ of integration around the branch point to = loq is equal to — ojq/A = — M /(4y / a) 
(see Appendix A). It should be noted that Ii is exactly equal to the first term in formula ( |2.3U| ) 
The sum of integral I2 along the interval (— iA, iA) of the imaginary axis and I± is 



for E, 



(2) reg 
C 



E, 



(2) reg. 



A 



c (R,A) = -—\n 



sinh RJA 2 + ul 



— dy\n 
2nJ y 



sinh RJ y 2 + lUq 



^. (3.15) 



Integration by parts is already done here. Formula 
00, A). From ( |3.15| ) it follows that 



4[) requires an asymptotics E^ reg (R 



E, 



(2) reg 
C 



(R -> 00, A) 



A 
2^ 



In 



sinh R J A 2 + ul 



a 

— jdy \Ry/tf + J*-hi2 



4 ' 
(3.16) 

The constant term — uJq/A is preserved here to satisfy condition (|3.1| ) which defines the behavior 
of the string potential at infinity. Otherwise this term would appear in the right-hand side 
of (|3.1|) , but that is physically unacceptable. At large distances string potential should be 
determined by its classical value only. Inserting ( |3.15| ) and ( |3.16| ) into ( |3.4| ) we find for i = 2 



E. 



(2) ren 
C 



(R) 



2tt 



du In ( 1 



R 



71 



u 2 du) 



^/uj 2 + u 2 e 2R V^!> 



(3.17) 



It is interesting to the compare the formula derived with an analogous expression for 
E ( c )ren (R) (see Eq. (3.9)). Formula (3.17) can be obtained by changing the variable fre- 
quency in (3.9) to \Juo 2 + uJq. This completely corresponds to the fact that Eq. ( |3.9| ) deals with 
oscillations of the massless (two-dimensional) scalar field on the segment [0, R] while Eq. (|3.17| ) 
treats oscillations of the same field, but with the mass equal to uq = M / yfa (see field equations 



At first sight, the expression obtained for Eq ren (R) by making use of the consistent renor- 
malization of the string tension does not coincide with that derived by analytic continuation 
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of the Epstein- Hurwitz zeta function (see formula fl2.31| )). This is not true, however. Equa- 
tions ( |3.17|) and ( |2.31| ) are completely equivalent. To show this, let us expand the logarithm 
in® 

Eg ren (R) = -— Y^n- 1 / due-W^. (3. 18) 

271 " =1 

By changing the variable u = u sinh t, the integral is reduced to the tabular one [16] 

00 

K v {z) = J e- zcosht cosh(ut)dt 


with z = 2nuoR. Finally we deduce the series Q2.31|) from (|3.17|) 



4 2) ren (R) = ^ fduln (l - e- 2 V^f) = £ n^K^nuoR), (3.19) 
n=1 

where Uq = M /y/a. Thus we found an integral representation for the series (|2.31|) . This series 
is convenient for investigating the behavior of the Casimir energy Eq ren (R) at large distances. 
Taking into consideration ( |2.28| ) we get 

E^ ren (R) ~_I^V /2 e -2-o/?. (3.20) 



R->oo 4\ttR 



The integral representation (3.17) enables one to study the asymptotics of Eq ren (R) at small 
R. From |TT7|) it follows that E^ ren (R) has a singularity when R = 0. For small R the main 



contribution to this integral is given by large to, therefore one can neglect here the dependence 
on u)q. This immediately gives the asymptotics of E^ ren (R) for R — > 



j^(2) ren ( 



00 

(R) ~ — / du In fl - e 2ujR ) = (3.21) 





Thus, consistent regularization of the divergent series (|2.12[) and subsequent renormalization 
of the string tension justify the rejection of the singular (pole) term and i?-independent constant 
in Eq. (|2.30| ) when analytic continuation of Epstein-Hurvitz ^-function is used. It is worthwhile 
to emphasize an important advantage of the proposed regularization by contour integration and 
subsequent subtraction. In this way we obtain the spectral representation for string energy at 
zero temperature (see Eqs. ( p.9|) and ( p. 171 ) in contrast to analytic continuation of ^-functions 
(Eqs. ( [2.20|) and Q2.31|) ). Proceeding from this spectral representation one can immediately 
derive the string free energy at finite temperature. To this end one must pass from integration 
to summation over the Matsubara frequencies fl n = 2imT, n — 0, ±1, ±2, . . .. Practically it is 
done by the substitution 

00 

duo -»• 2TxTduo^^ ~ ^™), (3.22) 

n=0 

where T is the temperature (see Appendix B). The prime of the sum sign means that the term 
with n = should be multiplied by 1/2. 
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For either quantity E^)(R) i = 1,2 we have obtained two integral representations (see 
Eqs. ( |3.9|) and (|3.17| )). Substitution (3.22) in these formulas with logarithmic functions gives 
us the free energy at finite temperature (Appendix B). For example, 



F^{R,T) = 2T£'m (l - e" 2 ^ . (3.23) 

n=0 ^ ' 

Taking the limit u — > in (3.23) one can obtain the free energy F^'(R, T) that diverges due to 
the term with n = (see Appendix B). Making the substitution (3.22) in the second version of 
the spectral representations (3.9) and (3.17) we arrive at the internal energy at temperature T 

oo 

UU(R,T) = -AnRT'Y:' tt " 7 , (3-24) 

^ exp(AirnRT) - 1 

U®(R,T) = -87r 2 RT 3 Y'- U ; . (3.25) 

n=o Jni + ^ 2 exp(2i? v /fi2 + W 2) _ j 



Both the energies, U^(R,T), i — 1,2 are well defined. The last two equations prove to be 
convenient for investigating the behaviour of the internal energies at large and small T. Let us 
demonstrate this using Eq. (3.24). At large T the main contribution to (3.24) comes from the 
first term with n = 

U {1) {R,T^ 00) = --. (3.26) 
At small T the Euler-Maclaurin formula 

00 op 

E'/H = J f(x)dx--f'(0) (3.27) 

n=0 

can be used. In the case under consideration 



As a result, we obtain for small T 



^Wr>— 5^-^. (3.29) 



4 Conclusion 

The experience of treating the divergences shows that a correct result can be obtained by 
applying practically any regularization and renormalization procedures provided that the pre- 
scriptions are properly modified. Therefore, when evaluating such methods, those should be 
preferred which are closer to the quantum field theory. Only in the framework of this approach 
one succeeds in formulation of a consistent renormalization procedure. Besides, quantum field 
formalism provides a clear and simple transition from zero temperature calculations to those 
at finite temperature [23]. In view of this, contour integration has an obvious advantage. At 
first it was proposed as a simple method for calculating the van der Waals forces between di- 
electrics [24] (see also [13, 17, 25, 26]). However its relation to the formalism of the Green's 
functions is not elucidated properly. And this problem is undoubtedly worth investigating. 
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Appendix A. Investigation of the contour integral deter- 
mining Eq re9 (R, A) 

Let us consider integral (3.14) for individual parts of contour C shown in Fig. 3. Integration 
along the semicircle of radius A contributes only to the counterterm, therefore we do not analyze 
it here. When going from the upper edge of the cut to the lower one, the integrand does not 
change. As a result, integrals along these two parts of the contour C are cancelled mutually due 
to opposite directions of integration. Only integration around the branch point and along the 
interval of the imaginary axis (— iA, iA) lead to finite contributions. While integrating around 
the branch point uj = ujq we introduce usual variables uj — ujo = pe 1 ^ with p — > 0, and in terms of 



them we have uj 2 — uJq = (uj + uj q )(uj — uj ) ~ 2uope l< ^, cos(RJuj 2 — uj" 2 ) ~ 1, sm(RJuj 2 — ujq 



Ryuj 2 — uj 2 .. Taking this into account we deduce 

2tt 



R 7 Jjpe^i d<p _ u _ M 
1 Am J 2u pe i <t>R 4 Ay/a' 1 ' ' 



o 



The integral I\ is exactly equal to the first term in ( |2.30| ) which is independent of R. When 
integrating along the imaginary axis, trigonometrical functions in (3.14) become hyperbolic 
ones 



R_ h-y 2 )idy coshCRyV + ^o 2 ) = A 
4 ™ / iJy 2 + ujI smh(Rjy 2 + ojI) V ^ { 



In sinh f Ry y 2 + Uq 

Summing (A.l) and (A. 2) and integrating by parts one arrives at formula (3.15) 



(A.2) 



Appendix B. Transition to finite temperature in infinite 
system of noninteracting oscillators 

Let us consider an infinite system of noninteracting oscillators with eigenfrequencies u n , n = 
1,2,... determined by the equation 

f(u,R) = 0. (B.l) 

Roots of this equation are assumed to be situated on the real axis in the complex plane uj. 
This set of oscillators arises, for example, in quantization of a scalar field defined on the line 
segment [0, R\. Boundary conditions imposed on this field result in frequency equation (B.l). 
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Without loss of generality, for relatistically invariant system one can admit that the function 
/ satisfies the condition 

f(-u,R) = f(u,R). (B.2) 
The free energy of this system is given by 



r 

F(R,T) = J2 ^+T\n(l-e 



n=l 



-ui„/T 



(B.3) 



At zero temperature this formula obviously turns into the energy of zero point oscillations 

In the general case sum ( |B.3| ) diverges, therefore to obtain the finite value for the free energy, 
we have to use the renormalization procedure discussed in Section 3. First, the infinite sum 
should be represented as the contour integral 



n=l 



1 / f J'(z,R) 
2m J f{z,R) 



(B.4) 



where 



UJ 



+ rin(l-e" 



■u/T 



CO 

T 



TJ2 

n=l 



-nuj/T 



n 



The contour C, as in Section 3, surrounds the first N roots of Eq. ( |B.1| ). 

As shown in Section 3, only integration along the imaginary axis gives a finite contribution 
to the free energy 



F re9 (R,A) 



1 



2m 



W 
2 



CO 1 \ 

TY,-e- my/T ) d\n[f(iy,R)}. 



n=l 



n 



(B.5) 



On integrating by parts we obtain 



F re °(R, A) = i jdy (i + f) e-^'A \n[f(iy, R)] 

_^ \ n=l I 



(B.6) 



-A *■ n - 1 



cos 



T 



i sin 



T 



The off-integral terms are omitted in (|B.6|) because they contribute only to the counterterm. 
Taking into account (|B.2|) we can drop terms with sine functions 



1 f 

F re °(R,A) = - dy 

7T J 



1 ^ (ny 

z n=l v 1 



Hf(w,R)} 



(B.7) 



The renormalized free energy is obtained by the subtraction 

F ren {R,T) = [F reg {R, A) - F re9 {R^ oo,A)}\ A ^ oc = 

f{iy,R) 



dy 



1 f (ny 

z n=l v 1 



hi 



f(iy, oo) 



(B.8) 
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With allowance for the Fourier-series representation of the ^-function 



]T 5(y -2nnT) = - + cos (Jjr 

n=—oo n=l 



integration in 



can be done to produce 

~ f(2mnT, R) 



F ren (R,T) =T ln 



f(2irinT, oo) 



n=0 



f(2mnT, R) 
f(2irinT, oo) 



(B.9) 



Now we apply formula ( |B.9|) to the models considered in Section 3. In the case of a scalar 
field with mass uq on the segment [0, R] we have 



and Eq. (|B.9|) gives 



f(u, R) = sin RJuj 2 + Wq 



FW(R,T) = 2T^'ln (l - e -*V<M) 

n=0 ^ 



(B.10) 



where Q n = 2imT. The same result was obtained in Section 3 by transition from the integral 
representation for the Casimir energy at zero temperature to summation over the Matsubara 
frequencies (see Eq. (3.23)). Proceeding from flB.lOQ one can derive the internal energy of the 
system under consideration applying thermodynamic rules 



U (2 \R,T) = — T 2 



d F®(R,T) 



dT 



T 



-8n 2 RT 3 J2' 



n 



o Jill + ul exp(2RjQl + ul 



This equation was derived in Section 3 by a simple substitution (see Eq. (3.25)). 

In the case of massless scalar field (ujq — > 0) the term with n = in QB.10|) diverges 



(B.11) 



F^(R,T) =2Tj2'\n(l-e 



n=0 



-AnnRT 



T lim 



exp(—A7mkRT) 



k=l 



oo 

2T In (l - e 



AnnRT 



n=l 



OO 1 oo 

= -T T + 2T In (l - e~ 4wnRT ) . (B.12) 

k=l * n=l 

This divergence is a manifestation of the well-known infrared instability of a massless scalar 
field in two-dimensional space-time. In Section 2 some reasons were given to attribute the 
Euler constant value, 7 to the sum of the divergent series Y^kLi k~ l ■ Finally we obtain a finite 
expression for the free energy of the massless scalar field on the segment [0, R] 



00 

F (1) (R, T) = — 7T + 2T In (l - e 



AwnRT 



n=l 



(B.13) 
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It should be noted that this treatment of infrared divergences in the problem in question is 
absolutely formal and it needs the physical justification. 

However the internal energy of this field is well defined. Putting u = in (B.ll) we get 

oo 

UW(R,T) = -AtiRT 2 ^' ™ . (B.14) 

In Section 3 the same formula has been derived by a formal substitution (see Eq. (3.24)). 
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Figure Captions 

Fig. 1. Contour C used in analytic continuation of the Riemann zeta- function. 
Fig. 2. Transformation of the contour in integral (3.6). 
Fig. 3. Contour used for summing the roots of Eq. (2.9). 
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